Feshbach Spectroscopy of a Shape Resonance 
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We present a new spectroscopy technique for studying cold-collision properties. The technique is 
based on the association and dissociation of ultracold molecules using a magnetically tunable Fesh- 
bach resonance. The energy and lifetime of a shape resonance are determined from a measurement 
of the dissociation rate. Additional spectroscopic information is obtained from the observation of a 
spatial interference pattern between an outgoing s wave and d wave. The experimental data agree 
well with the results from a new model, in which the dissociation process is connected to a scattering 
gedanken experiment, which is analyzed using a coupled-channels calculation. 
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The field of ultracold molecules has seen tremendous 
progress in the past two years. The key technique that 
triggered this development is the association of ultracold 
molecules from ultracold atoms using a Feshbach reso- 
nance HEU&HEHHIi. A Feshbach resonance is 
caused by the resonant coupling of a colliding atom pair 
to a molecular bound state. The system can be tuned into 
resonance by applying a magnetic field. To produce sta- 
ble molecules, the magnetic field is slowly ramped across 
the Feshbach resonance to lower fields. For detection, the 
molecules are dissociated by ramping the magnetic field 
back to higher values above the Feshbach resonance. In 
this field range, the energy of the molecular state lies 
in the continuum of unbound atom pair states. The 
coupling between the discrete molecular state and the 
continuum of atomic states leads to an exponential de- 
cay of the molecules into unbound atom pairs. These 
atoms are then imaged with standard techniques. Be- 
yond this molecule detection scheme, the dissociation has 
other interesting applications: For example, the width of 
the Feshbach resonance can be extracted from the energy 
released in the dissociation process 0, El ■ 

In this letter, we use the dissociation process as a new 
spectroscopic tool for studying cold-collision properties. 
Specifically, we investigate a d-wave shape resonance in 
87 Rb. The shape-resonance state is a quasi-bound state 
which is localized behind the centrifugal barrier. The 
energy of this state lies in the continuum of unbound 
atom pair states. In our experiment, the shape reso- 
nance is probed by tuning the energy of the Feshbach 
molecules with the magnetic field. If the energy of the 
molecules matches the energy of the shape resonance, 
then population is transferred from the molecular state 
to the shape-resonance state. From there, the popula- 
tion decays rapidly into unbound atom pairs by tunnel- 
ing through the centrifugal barrier. The molecules can 
therefore dissociate in two ways, either directly into the 
continuum or indirectly by passing through the shape- 
resonance state. The indirect processes turn out to be 
much faster than the direct ones, so that the molecule 
dissociation rate is drastically enhanced. However, if the 



energies of the molecules and the shape resonance do not 
match, then indirect processes will be suppressed. The 
magnetic-field dependence of the dissociation rate thus 
reveals the energy and lifetime of the shape resonance. 
This information can be used as spectroscopic input to 
constrain theoretical models for the cold-collision prop- 
erties. Little is gained in the case of 87 Rb, where the 
collisional properties are rather well known, and where 
the shape resonance has been observed in previous ex- 
periments [ill H2I fl3| . However, our method is inter- 
esting for systems, where the collisional properties are 
unknown, such as heteronuclear mixtures or non-alkali 
species. Here, theoretical models starting from ab-initio 
calculations presently do not have enough experimental 
input to obtain realistic predictions for the cold-collision 
properties. 

The ex per imental setup was described in detail else- 
where H, llfl fl4| . In brief, a Bose-Einstein condensate 
(BEC) of 87 Rb atoms in the hyperfine state \f, mt) = 
|1, 1) is prepared in an optical dipole trap. The magnetic 
field is held slightly above the Feshbach resonance, which 
is located at B res ~ 632 G 0] . It is caused by a d-wave 
molecular state and has a width of 1.3 mG The 
BEC is released from the trap and 1 ms later molecules 
are created by slowly ramping the magnetic field down- 
ward across the Feshbach resonance. Remaining atoms 
are spatially separated from the molecules using a Stern- 
Gerlach field. Next, the molecules are dissociated back 
into unbound atom pairs by jumping the magnetic field 
to a value B > B res and holding it there for a variable 
time thoid 0| • During thoid, population in the molecu- 
lar state decays exponentially with a rate r(B), which 
depends on the value of B during thoid- After thoid, the 
magnetic field is switched off rapidly, which stops the 
dissociation process. After 0.85 ms time of flight (count- 
ing from the dissociation), an absorption image is taken. 
Molecules that did not decay during thoid are invisible 
in the image, because only unbound atoms resonantly 
absorb light from the detection laser beam. 

The dissociation rate T(B) is extracted from a series of 
images for fixed B and variable thoid- The number of dis- 
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FIG. 1: Molecule dissociation rate F as a function of magnetic 
field B with respect to the position of the Feshbach resonance 
B res . The experimental data for the total rate (circles) clearly 
show the effect of the shape resonance near 1.3 G. The solid 
line is a Lorentzian fit to the data. The dashed (dotted) line 
shows the prediction for the total (s-wave) dissociation rate 
from a coupled-channels calculation. The s-wave dissociation 
is not affected by the shape resonance. 



sociatcd atoms is determined from each image and then 
fit to an exponential as a function of thoid- This yields 
the dissociation rate shown in Fig. ^ The enhanced dis- 
sociation rate near 1.3 G due to the shape resonance is 
clearly visible. 

In a simple analytic model, the dissociation can be 
described as a two-step process. In the first step, the 
interaction Hamiltonian H transfers population from 
the molecular state \ip mo i) to the shape-resonance state 
l^shape)- This is described by a generalized Rabi fre- 
quency £1 = (2/h)\(tp shape \H\'ip mo i)\. In the second step, 
the population tunnels with a rate T shape from \tp shape) 
into the continuum of unbound atom pair states. One 
can show that for fl <C r s /p. the molecular state decays 
exponentially with a rate |l6| 
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where E s hape and E mo i are the energies of the shape- 
resonance state and the molecular state, respectively. 
E mo i depends on B, while T shape , E shape , and ft in this 
simple model do not. For the small magnetic- field range 
considered here, a linear approximation holds 



E mo i{B) = (B — B res )Afx , 



(2) 



where B res is the position of the Feshbach resonance 
and A/z is the difference in the magnetic moments be- 
tween a molecule and an unbound atom pair. A/i is de- 
termined experimentally from the time-of-flight images 
shown in Fig. [2J The radial cloud size in the images 
yields the kinetic energy E in the relative motion of the 
dissociated atom pairs. Due to energy conservation in 
the dissociation process, E = E mo i(B). We thus obtain 
Afi = k B x 230(18) zzK/G. Alternatively, Afi can be de- 
termined with a Stern-Gerlach method [3J. 
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FIG. 2: (Color online) Time-of-flight images of unbound 
atoms obtained in the molecule dissociation. The images were 
taken at different values of B — B res . The magnetic field is 
vertical in the images. The interference between the s and d 
partial wave undergoes a change in relative phase and ampli- 
tude. At 0.1 G, the dissociation is mostly s wave, producing a 
circle. For higher magnetic fields, both partial waves are pop- 
ulated. From 0.2 to 0.6 G, the interference suppresses atom 
emission along B, whereas the opposite relative phase in the 
interference between 1.0 and 1.4 G suppresses emission per- 
pendicular to B. At 0.7 and 0.9 G, the relative phase is such 
that neither component is strongly suppressed. The typical 
radius reached by the atoms during the constant time of flight 
increases for increasing B, thus indicating an increase of the 
kinetic energy released in the dissociation. 



T(B) in Eq. QJ is a Lorentzian, which is fit to the 
data as shown in Fig. 2] Using the above value for A/z, 
the best-fit parameters are E s hape — X 312(25) /zK, 
T sh ape = 15(3) MHz, and ft = 2ttx0.61(7) MHz. We thus 
extract the energy and lifetime of the shape resonance as 
well as the Rabi frequency without making any reference 
to the coupled-channels calculation further below. 

Figure [21 shows time-of-flight images of the atoms after 
dissociation. They exhibit spatial interference patterns 
created by different outgoing partial waves. Only two 
partial waves contribute to these interference patterns, 
namely those with quantum numbers (I, mi) = (0,0) and 
(2, 0) for the rotation of the atoms around each other 
|17|. The indirect dissociation processes (where popu- 
lation passes through the shape-resonance state) create 
outgoing d waves, because the shape-resonance state is 
a d-w&ve state and the tunneling out of this state does 
not affect the angular momentum. In contrast, the di- 
rect dissociation processes are independent of the shape- 
resonance and preferentially populate the s-wave state, 
because all other partial waves are suppressed for low 
energy by the centrifugal barrier. Note that interference 
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between different partial waves of cold atoms has_ 
ously been observed in scattering experiments |~ 
Previous molecule dissociation experiments, however, ob- 
served only outgoing s waves. 

The interference pattern is described by the following 
wave function 

^decayif) = 9 (r, t) (^A^OO - ^ yfihYmtf)) (3) 

where Yi ra (§ 1 ip) is a spherical harmonic, the real num- 
bers Pi and 5 re i characterize the amplitudes of the par- 
tial waves and their relative phase, and g{r, t) is a radial 
wave function. g(r, t) is normalized such that P0+P2 = 1, 
which implies that /?; is the branching ratio for decay into 
the l-th partial wave. 

The branching ratio and the relative phase S re i are 
extracted from the images with the same method as used 
in Ref. [l3j . In brief, computed tomography is used to 
reconstruct the 3-dimensional (3D) density distribution 
from the 2D images. The 3D density is sorted into 20 
bins to obtain the probability W(i3) for finding an atom 
at angle A fit of the angular part of the modulus 
squared of Eq. J2J to W(t?) yields the fit parameters #2 
and S re i . These parameters are shown in Fig. |21 The 
branching ratio clearly shows the enhanced decay into 
the d wave due to the shape resonance. For (3^ close to 
or 1, the fit cannot reliably determine the relative phase. 

In order to analyze the dissociation more rigorously, 
we performed a coupled-channels calculation for a scat- 
tering gedanken experiment that is closely related to our 
dissociation experiment. We briefly summarize the the- 
ory here. For more details see Ref. |l9). In the gedanken 
experiment, we consider a colliding atom pair with ki- 
netic energy E in the relative motion in the presence of 
a magnetic field B near the 632 G Feshbach resonance. 
(S-matrix elements Su> are calculated in the _E-_B-plane 
for I and V equal to or 2 (both with mi = 0). Accord- 
ing to the theory of multichannel scattering resonances, 
these ^-matrix elements are ex pec ted to follow an ana- 
lytic expression of the form [jfl 12 lj 
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Here, 6w is the Kronecker symbol, T(E) is the total decay 
rate of the molecular state, E mo i(B) is given by Eq. J2J, 
and the significance of the partial-decay-rate parameters 
Tu'(E) will be discussed below. Finally, S^ 9 (E) is the 
background scattering phase for the l-th. partial wave, 
where background means for magnetic fields far away 
from the Feshbach resonance 

For any given E, we fit Eq. (@J for variable B to the 
5-matrix elements obtained from our coupled-channels 
calculation. The fits match the ^-matrix extremely 
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FIG. 3: Parameters extracted from the spatial interference 
patterns, (a) Branching ratio for decay into the d wave. 
Again, the d-wave shape resonance is clearly visible, (b) 
Partial- wave phases and relative phase. The solid, dashed, 
and dotted lines show the theoretical prediction for the rel- 
ative phase S re i = S^ 9 — 8q 9 , the s-wave phase S^ 9 , and the 
d-wave phase S^ 9 , respectively. Experimental data (circles) 
agree well with the theory (solid line) in (a) and (b). 



well. The open-channel physics, including the shape res- 
onance, is independent of B and is therefore included in 
the energy-dependent fit parameters 5^ 9 (E) and Tui(E). 
Furthermore, the combination of the fits yields A/i = 
k B x 224 iM/G. 

We will now discuss the connection between our dis- 
sociation experiment and the fit parameters S^ 9 (E) and 
Tu>(E). As mentioned earlier, the kinetic energy E of the 
atom pairs after dissociation is E = E mo i(B). Therefore 
the energy-dependent fit-parameters for the gedanken ex- 
periment become magnetic-field dependent parameters in 
the dissociation experiment 5\ 9 (B) and Tw (B) . 

The link between the scattering gedanken experiment 
and the dissociation experiment can be established from 
evaluating the different terms in the asymptotic form of 
the regular scattering wave function 
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where we assume in the gedanken experiment that only 
one partial wave V is initially populated. Here, k is re- 
lated to the energy by E — h 2 k 2 /(2m re d) with the re- 
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duced mass m re d, and we split the 5-matrix into a back- 
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ground part and a resonant part Sw 

S h l § = e 2ish i 9 W5 lv . 

The scattering wave function consists of three parts: 
an incoming wave, an outgoing background wave, and 
an outgoing resonant wave. The resonant part is due 
to particles that make the transition to the molecular 
state and subsequently decay back into the open channel. 
The Stern-Gerlach separation removes all incoming flux 
and along with it the background scattered wave. Hence, 
these two terms need to be removed from the scattering 
state in order to describe the decay wave function 
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Using the above equations, this can be rewritten to 
yield Eq. where the relative phase of the decay- 
ing partial waves is 5 re i 
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ing ratio for decay into the Z-th partial wave is given by 
= T U {B)/T{B). This means that T U {B) repre- 
sents the partial decay rate into the /-th partial wave. 
The theory for the decay rate, branching ratio, and the 
phases are shown as lines in Figs. ^ an d El The good 
agreement between theory and experiment without any 
free fit parameters supports the theory developed here. 

Finally, we extract the parameters of the simple model 
from the full theory. The d-wave decay rate = T — 
Too is the difference between the total and the s-wave 
decay rate, both shown in Fig.^ The maximum of is 
located at 1.28 G corresponding to E s h ape = /sbx287^K. 
The value and the curvature of (B) at the maximum 
correspond to r s h Qpe = 17 MHz and Q — 2tt x 0.61 MHz. 
These numbers agree well with the values obtained by 
fitting the simple model Eq. Q to the experimental data. 

In conclusion, we showed that dissociation of ultracold 
molecules can be used to measure the energy and life- 
time of a shape resonance. Additional information was 
obtained from the analysis of the spatial interference pat- 
terns. This is enough information to constrain theoreti- 
cal models so much that the basic scattering properties 
of the system can be determined. We developed a new 
model for the dissociation, which agrees well with the 
experimental results. 
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